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Abstract 

With a view towards models of quantum computation and/or the inter- 
pretation of linear logic, we define a functional language where all functions 
are linear operators by construction. A small step operational semantic (and 
hence an interpreter/simulator) is provided for this language in the form of a 
term rewrite system. The linear-algebraic A-calculus hereby constructed is 
linear in a different (yet related) sense to that, say, of the linear A-calculus. 
These various notions of linearity are discussed in the context of quantum 
programming languages. 

1 Introduction 

Quantum computation lacks a convenient model of computation. To this day its 
algorithms are expressed in terms of quantum circuits, but their descriptions al- 
ways seem astonishingly remote from the task they do accomplish 1 12 1. Moreover 
universality is only provided via the notion of uniform family of circuits |28|. 
Quantum Turing machines solve this latter point, yet they are even less suitable 
as a programming language 0. Another approach is to enclose quantum cir- 
cuits within a classical imperative-style control structure 112.01 — but we wish to 
avoid this duality, in an attempt to bring programs closer to their specifications. 
Functional-style control structure, on the other hand, seem to merge with quan- 
tum evolution descriptions in a unifying manner. With a view towards models 
of quantum computation, we describe a functional language for expressing linear 
operators, and linear operators only. 

We are careful, however, not to bury our presentation of this language of linear 
operators within too many quantum computation-specific considerations. The aim 
is to reach an audience of logicians also, as we suspect a strong connection with 
issues of computational interpretations of linear logic. 
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We provide a semantic for the language in the form of a term rewrite system 
II II . These consist in a finite set of rules I — ► r, each interpreted as follows: 
"Any term t containing a subterm al in position p (i.e. t — t[al] p ) should be 
rewritten into a term t' containing or in position p, with all the rest unchanged 
(i.e. t' = t[ar] p )". Here a denotes a variable substitution. The minimalist in- 
terpretation of the rules makes term rewrite systems (TRS) extremely suitable for 
describing the behavior of a computer languages unambiguously — so long as 
the order in which the reductions occur does not matter to the end result (a prop- 
erty named confluence). Moreover, because I — ► r may be seen as an oriented 
version of equation I = r, the TRS provides both an operational semantic (an 
interpreter/simulator for the language) and an axiomatic semantic (an equational 
theory in which to prove properties about the language). 

We begin with a simple language for vectors containing constants for base 
vectors, addition and product by a scalar. On terms of this language we define 
a rewrite system reducing any term expressing a vector to a linear combination 
of base vectors [4|. We have also proposed in |4| an extension to a language 
containing a tensorial product operation (Section|2j. 

Such a language must rely on a language and rewrite system for scalars. This 
raises the problem of the conditional rewriting required for division, which we can 
circumvent, basing quantum computation upon the ring of diadic floats together 
with -j= and imaginary number i (Section [3}. More generally, it should be said 
that a language of linear operators does not need division. 

Modern days functional languages such as Caml, Haskell etc. are based upon 
two basic evaluation mechanisms: matching, which provides conditional branch- 
ing by inspection of values; and some avatar of the X-calculus. The first mecha- 
nism is obtained as we extend the term rewrite system to handle linear maps — 
themselves denoted as superpositions of bipartite states, e.g. 

(true > false + false d> true) * false — >* true. 

Applications are therefore analogous to contractions in tensorial calculus: this 
approach offers an elegant paradigm to represent quantum operations as quantum 
states (Section^. 

The second mechanism is obtained through an implementation of A-terms via de 
Bruijn indices, a scheme whereby variables are encoded as integers referring to 
their binders, e.g. 

\x.(Xy.(x <g> y)) is encoded as L(£(var(l) ® var(O))). 
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The question of the interpretation of terms such as Xx.(x (g> x) is lengthily ad- 
dressed as we draw a distinction between cloning and copying. The semantic of 
our calculus forbids only the former, non-linear operation, by enforcing a higher 
priority of the addition's distributivity over substitution (Section 0. This is fol- 
lowed by a short example taken from our implementation (Section|6j. 

Erasure on the other hand remains allowed in our calculus, because we do 
not restrict ourselves to unitary operations. Whilst we discuss possible well- 
formedness conditions to implement this restriction (a crucial one for quantum 
computation), the claim here is to have provided a "linear" A-calculus, in the 
sense of linear algebra. We discuss the various notions of "linearity" used in 
quantum programming languages, such as the one by Van Tonder |23| (Section 

m 

2 Vectorial spaces 

We seek to represent quantum programs, their input vectors, their output vectors 
and their applications as terms of a first-order language. Moreover we seek to 
provide rules such that the term formed by the application of a quantum program 
onto its input vector should reduce to its output vector. Several terms may be 
used to express one output vector, as a consequence we must ensure that these all 
reduce to one unique, normal form, upon which there is nothing more to compute. 
The most natural normal form to aim for is that of a linear combination of the base 
vectors, i.e. the computation finishes once we have the coordinates of the output 
vector. 

We start with the language of vectorial spaces, i.e. a two-sorted language £ 
having sort K for scalars and sort E for vectors — together with: two constants 
and 1 of sort K; a constant of sort E; two binary symbols + and x of rank 
(K, K, K); a binary symbol + (also) of rank (E, E, E); and a binary symbol . 
of rank (K, E,E). In |4| we described a term rewrite system reducing any term 
expressing a vector into a linear combination of base vectors. The term rewrite 
system develops 

4. (false + true) — ► 4.false + 4. true 

but factorizes 

4.false + 6.false — > (4 + 6). false. 

according to the rules in figure^ Such a TRS arises as we orient six of the eight 
equations axiomatizing vectorial spaces. Only those two axioms corresponding to 
associativity and commutativity of vector addition are left aside, because we use 
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rewriting modulo AC(+). Moreover we need to add three more rules for conflu- 
ence. 



Figure 1 : Vectorial spaces 


A.(u + v) — ► 


A.u + A.v 


A.u + /i.u — ► 


(A + /i).u 


X.(fi.u) — ► 


(A x /i).u 


u + — ► 


u 


l.u — > 


u 


O.u — > 





A.u + u — ► 


(A + l).u 


u + u > 


(l + l).u 


A.O — ► 





with + an AC symbol. 



But these rewrite rules do not take into account computation on scalars. The 
latter must be added by mixing in another rewrite system S, rewriting scalar to a 
normal form. 

Definition 1 (Scalar rewrite system) A scalar rewrite system is a rewrite system 
on a language containing at least the symbols +, x, and 1, such that: 

• S is terminating and ground confluent, 

• for all closed terms A, p, and v, the pair of terms 

- + A and A, 

- x A and 0, 

- 1 x A and A, 

- A x (jit + v) and (A x /z) + (A x v), 

- (A + /u) + v and A + (fi + v), 

- A + n and fi + A, 

- (A x fj,) x v and Ax(/ix v), 

- A x /i and fi x A 
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have the same normal forms, 

• and 1 are normal terms. 

The following propositions can be found in @. 

Proposition 1 For any scalar rewrite system S, the rewrite system RU S is ter- 
minating and ground confluent. 

Proposition 2 If t is a normal close term whose constants are amongst Xi , . . . , x„. 
The term t is or it has the form \\Xi 1 + ...Xk^i k +Xi fc+1 +x JJc+i where the indices 
%i, ik+i are distinct and the A& 's are neither nor 1. 

Note that the algorithm defined by R is relatively common in computing, for 
presenting any vector as a linear combination of base vectors. But it does in fact 
define vectorial spaces, as any mathematical structure validating the algorithm. In 
this sense we have provided a computational definition of vectorial spaces. 

Furthermore note that the support for tensor products is easily added into the 
TRS, through the six rules given in figure|2] PropositionlOremains true when R is 
extended with those six additional rules, whilst proposition|2]now yields normal 
forms for terms in E <£> E of the form or 



AiXi! <g>y n + ... + A fe x ifc (g> y jk - 


1- x lfc+1 <g> y ]k+1 + ...+ x lfe+ , <g> y Jk+l , 


where the pairs of indices («i , Ji) , ■ ■ • 


, (ik+i,jk+l) are distinct and the A fc 's are 


neither nor 1 1 4 1 . 




Figure 2: Vectorial spaces: tensors 


(u + v) ® w - 


— >U®W + V(g)W 


(A.u) ® v - 


— > A.(u(g) v) 


u <g> (v + w) - 


— >U(g)V + U(g)W 


u® (A.v) - 


— > A.(u(g) v) 


0(g) u - 


— > 


u® - 


— » 
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3 The field of quantum computing 



Fields are not easily implemented as term rewrite systems, because of the con- 
ditional rewriting required for the division by zero. In the previous section such 
problems were avoided by simply assuming a TRS for scalars having a certain 
number of properties, but if the objective is to lay the ground for formal quantum 
programming languages, then we must provide such a TRS. The present section 
briefly outlines how this is achieved. 

3.1 Background 

We seek to model quantum computation as a formal rewrite system upon a finite 
set of symbols. Since the complex numbers are uncountable, we must therefore 
depart from using the whole of C as the field K of our vectorial space. Such 
considerations are commonplace in computation theory, and were successfully 
addressed with the provision of the first rigorous definition of a quantum Turing 
machine 1 6 1 . In short the quantum Turing machines are brought as an extension 
of probabilistic Turing machines 

IQ: head states, S: alphabet, 
5: transition function, q ,qf ■ start,end stat$ 

whose transition functions are no longer valued over the efficiently computable 
positive reals (probabilities) 

<5 : Q x £ — > (Q x £ x {Left, Right} -> M+) 

but over the efficiently computable complex numbers (amplitudes) 

<S : Q x £ — > (Q x £ x {Left, Right} -► C). 

In both cases S is constrained to be a unit function (probabilities/squared modulus 
summing to one), and for the quantum Turing machine S is additionally required 
to induce a unitary global evolution. A well-known result of complexity theory is 
that probabilistic Turing machines remain as powerful when the transition func- 
tion S is further restricted to take values in the set {0, |, 1}. The result in |6| 
is analogous: quantum Turing machines remain as powerful when the transition 
function 5 is further restricted to take values in the set {—1, — 0, -^=, 1}. Later 
it was shown in [2|, and independently in 1221 that no irrational number is nec- 
essary, i.e. 5 may be restricted to take values in the set {— 1, — ~, — |, 0, |, |, 1} 
without loss of power for the quantum Turing machine. 



6 



In the circuit model of quantum computation the emphasis was placed on the 
ability to approximate any unitary transform from a finite set of gates. This line 
of research (cf . 1 2 1 1 1 1 5 1 to cite a few) has so far culminated with [ 7 1 , where the 
following set 
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was proven to be universal in the above strict sense. A weaker requirement for a 
set of gates is the ability to simulate any unitary transform, a notion which is also 
referred to as encoded universality — since a computation on n qubits may for 
instance be represented as a computation on n + 1 "real bits", through a simple 
mapping. A recent paper shows that the gate 



G = 



/I 













1 










a 


-b 




V o 


b 


a 


) 


_ 3 
5 


and b 




8 
5' 



with either a = 6=^|, ora = | and b = |, has this property 1191 . Do 
appreciate how the result falls into line with those regarding the quantum Turing 
machine. 

Definition 2 We call computational scalars, and denote K the ring formed by the 



additive and multiplicative closure of the complex numbers {—1,1, 



V2' 



Once we have shown that the computational scalars arithmetics can be performed 
by a TRS, it will be sufficient to express the basic gates Q in our formalism to im- 
mediately obtain the more traditional notion of quantum computation universality. 
Hence our choice. 



3.2 Rules 

We begin by implementing natural numbers and unsigned binary numbers. That 
such TRS can be made ground confluent and terminating are now well-established 
results |9||27|. This places us in a position to build up diadic floats out of a sign, 
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Figure 3: Diadic floats 




fl(»,n::0,S(p)) 


— > fl(s,n,p) 




fl(neg,0,p) 


— ► fl(pos,0,p) 




fl(s,0,S(p)) 


— ► fl(s, 0,zeron) 


fl(pos, mi 


e\) timesf fl(neg, m 2 , e 2 ) 


— ► fl(neg, mi timesb m 2 , addn(ei, e 2 )) 


fl(neg,TOi 


ei) timesf fl(pos, m 2 , e 2 ) 


— ► fl(neg, mi timesb m 2 , addn(ei, e 2 )) 



an unsigned binary number and an exponent, e.g. fi(neg, 1, S'(zeron)) is to stand 
for — i, as exemplified in figure[3] 

Reached this point it suffices to notice that K, i.e. diadic floats together with 
imaginary number i and real number can be viewed as a four-dimensional 
module upon diadic floats. Indeed any such number could be represented as a 
linear combination of the form: 

a.l +(3.—= +7.i + <5.^=. 

As a consequence we can reuse the results of section |2] to implement computa- 
tional scalars and their additions. Computational scalars multiplication then needs 
to be defined, we do so modulo AC in figure 0] Notice that we overload the sym- 
bol x for multiplication of diadic floats and for multiplication of computational 
scalars. 

We conjecture that this TRS is ground confluent and terminating, but have not 
yet a formal proof for this assertion. 

Notice we have never defined a division operation. This is because only the 
ring properties of these numbers are required for expressing linear operations: we 
place ourselves upon a "module" rather that a full vectorial space. 

4 Matching construct 

We now turn to the definition of the matching constructs in our language. As we 
shall see, these constructs are nothing else than a reformulation of the rules for the 
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Figure 4: Scalar multiplication 

1 X v — ► V 

— x — — ► fl(pos, 1, 5(zeron)).l 

1 . i 

—= x i — > —= 

-j= x -j= — ► fl(pos, 1, S(zeron)).i 
i x i — ► fl(neg, l,zeron).l 
i x -j= — > fl(neg, l,zeron).-^ 

-j= x -i= — ► fl(neg, l,S(zeron)).l: 

(A.u) x v — > A.(u x v) 
(t + u) x v — > t x v + u x v 

with x an AC symbol. 



tensor product. 

4.1 Notations 

Your typical functional language (Haskell, ML. . . ) will always have "match- 
ing" constructs (for branching). For instance, here is a piece of Caml: 
let rec not b = match b with 
I false -> true 
I true -> false ; ; 

We wish to provide such constructs in our linear-algebraic calculus. Strangely 
enough these matching constructs are very close to the tensorial product con- 
structs. 

Mathematicians and physicist in this field would write linear maps instead: 
NOT = | true) (false | + | false) (true |. However here the (false | and (true| maybe 
viewed as patterns, waiting to be compared to the input vector through a scalar 



9 



product. Thus we choose to reconcile both worlds and write: 



NOT = false > true + true > false. 

An expression (t t> u) applied to a vector v will then reduce into (t • v).u, with • 
the scalar product. In this sense (t > u) * v does return u in so far as t overlaps 
with v. More formal justifications, and formal rewrite rules follow in the next two 
subsections. For now we give the reduction steps involved in the application of 
the phase gate P upon the vector true, as a motivating example for these rules: 

(false > false) + true > (-^= + i-^=).true^ * true 

— >* (false > false) * true + (true > (—= + i— r=).true) * true 

v 2 V 2 

— >* (false • true).false + (true • true).((— = + i— =).truel 

v V2 V2 ' 

, 1 1 X 

— ► O.false + (— = + i—j=)- true 

V2 y/2 

^* ( 7! + z 7I ) - true - 

All of the three gates forming a universal set for quantum computation are trivially 
expressed as terms in this notation: 

CNOT = (false <g> false) > (false ® false) 
+ (false (g> true) > (false ® true) 
+ (true ® false) t> (true ® true) 
+ (true ® true) t> (true <g> false) 

H = ^false > -^=. (false + true)^j + ^true d> -^=. (false — true) 
P = (false > false) + ^true t> (-^= + i-^=).true^ . 

4.2 Rules 

Since \> is just another type of tensor product, bilinearity applies (see figure [5] 
Notice the conjugation of the A scalar, denoted A, easily implemented in the TRS). 
Other than its left-hand-side antilinearity, the particularity of > is the reduction it 
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Figure 5: Matching operators bilinearity 



(t + u) > V 




t > V + U > V 


t > (v + w) 




t D> v + t D> w 


(A.u) > v 




A.(u > v) 


U > (fJb.v) 




(j,.(u > v) 


> u 







ut> 







(u + v) * w 




u * w + v * w 


(A.u) * v 




A.(u * v) 


u * (v + w) 




U * V + U * w 


u * (A.v) 




A.(u * v) 


* u 







u * 







(t + u) • V 




t • V + U • V 


t • (v + w) 




t • V + t • w 


(A.u) • v 




A.(u • v) 


u« (/x.v) 




/i.(u« v) 


• u 







u» 
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Figure 6: Matching operator and the scalar product 


(t t> u) * v — 


->(t.v) 


.u 


(t <g) u) • (v <g> w) — 


->(t.v) 


x (u • w) 


(t > u) • (v > w) — 


->(t.v) 


x (u • w) 


true • true — 


-> 1 




true • false — 







false • true — 


-> 




false • false — 


-> 1 





induces when placed left of an application symbol *, as described in figure[6] For 
definiteness, we may also add the rules in figureQ 

Notice that for now programming language does not use any variables. All 
functions are defined by adding elementary functions mapping base vectors to 
base vectors. And all functions are linear by constructions. Moreover until Section 
[?1 we do not worry about normalization and unitarity conditions. For now the the 
[> notation provides exactly what is needed: linear operations can be encoded as 
sums of tensor states describing which vector is associated to which 1 3 1 [ 8 1 . 

5 Lambda calculus construct 

We now turn to the last ingredient of our language: the lambda-calculus con- 
structs. In this case we will need to introduce bound variables. These bound 
variables will be handled through the use of de Bruijn indices. We shall also be 
more careful, as lambda-calculus usually allows to express non linear functions, 
whilst we have to restrict ourselves to linear ones. 

5.1 Cloning 

As we seek to provide a language of linear operator, we must pay attention to the 
fact that duplication of a vector can lead to nonlinear evolutions. For instance, 
whatever formalism we choose for quantum theory (vectors or density matrices), 
quantum operations act linearly upon their input states. This, in turn, implies that 
quantum states cannot be cloned. Indeed such an evolution acts upon a qubit as 
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Figure 7: Orthogonality rules 


(tgu)«(v>w) — 


-> 


(t <g> u) • true — 


-> 


(t <g> u) • false — 


-> 


(t [> u) • (v ® w) — 


-> 


(t [> u) • true — 


-> 


(t > u) • false — 


-> 


true • (v <g> w) — 





true • (v > w) — 


-> 


false • (v <& w) — 


-> 


false • (v t> w) — 


-> 



follows: 



(a|0)+/3|l»® 



(a|0)+/3|l»®(a|0)+j9|l)) 



/ 


a 






f ^ 


\ 









CLONE 


af3 






/? 






a/3 




V 





1 




V /3 2 


/ 





a 






/ a \ 









COPY 







/? 









V 





J 




V P I 



which cannot be linear (a more formal discussion can be found in |26|). Cloning 
should be distinguished from copying however, as we now illustrate once more on 
qubit: 

|0> <g> |0) ffi |0) ® |0) 

|i)®|o) c ^|i)®|i) 

(a|0)+/3|l)) ® |0) ffi a |o) g |0) +/3|1) ® |1) 
Such an evolution is perfectly valid, and in the above case it may be imple- 
mented as a single application of the quantum gate CNOT, which is of course 
both linear and unitary. 

In classical functional languages terms such as Xx.f(x, x), with (\x.f(x, x) t) — 
f(t, t), are crucial for the expressiveness. Recursion, for instance, relies upon such 
terms, and is absolutely necessary for universality. When designing a quantum 
functional language we therefore face a choice: 
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• Either we prevent terms such as Xx.f(x, x) from being applied to quantum 
states — thereby ensuring that no quantum cloning is allowed. But we must 
authorize their applications upon "classical terms" for expressiveness. As a 
consequence the language must be able to keep track of quantum resources 
versus classical resources. This is the approach followed by Van Tonder 



• Or we allow terms such as Xx.f(x, x) from being applied to quantum states 
— only to be interpreted as a quantum copy. We may still want to keep track 
of quantum resources versus classical resources, but not for the purpose of 
forbidding cloning. This is the approach we take. 

For now the latter option seems preferable, since it models the dos and don 'ts of 
the linearity requirement more closely, whilst keeping the calculus to a minimum. 
Moreover copying, as we now show, can be imposed over cloning by the semantics 
of the calculus alone. Thus, the fact that classical states can be cloned is proved 
and not postulated in our language. 

5.2 Substitution of de Bruijn indices 

The point of the previous discussion is that we can only duplicate basis vectors. 
Informally 

(\x.(x®x)) * true — >* true ® true is OK; 
{\x.{x®x)) * (false + true) — >* ((\x.(x®x)) * false) + ((\x.(xtgix)) * true) 



(\x.(x®x)) * (false + true) — >* (false + true) <g> (false + true) is not OK. 

Again another way to grasp this idea is to realize that faced with a term of the 
form (Xx.t) * (u + v), one could either start by proceeding to the substitution, or 
start by applying the right-hand-side linearity of *, leading to two different results. 
So that operations remain linear, we must favour the right-hand-side linearity of 
* over substitution. The rules of figure [8] accomplish exactly that. The three first 
rules are the most straightforward, they invoke the linearity of the vector to be 
substituted. The three following rules treat the base cases, when the vector to be 
substituted is down to a basic state. The last two rules handle the more subtle case 
of tensor states u ® v or u > v. In a word the trick is to treat 



|23|. 



(true <g> true) + (false <g> false) is OK; 
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Figure 8: Enforcing linearity over substitution 
L(u) * v — ► (u of v) 
t of (r.v) — ► r.(t of v) 
t of (v + w) — > (t of v) + (t ofw) 
t of true — ► t bof subst(true) 
t of false — > t bof subst(false) 
tof — ► tbofsubst(O) 
t of (v <g> w) — ► 

((((tbof fr(T))bof fr(t)) bofsubst(var(0)<g)var(S(0)))) ofv) ofw 
t of (v > w) — ► 

((((tbof ft (T)) bof lr(t)) bofsubst(var(0)t>var(S(0)))) ofv) ofw 



with y a fresh variable, and then proceed recursively. 

Once the vector to be substituted is a basic state, we can safely proceed to the 
substitution using a calculus of explicit substitutions ||T| 1101 1161 . Here we have 
chosen to represent avariables by their de Bruijn indices, i.e. each variable is now 
an integer number corresponding to number of binders ("L" or "A" symbols) one 
must go through before reaching the binding occurrence. For instance 

Xx.(Xy.(x ® y)) is encoded as L(L(var(l) ® var(O))) 

since there is one A symbol lying between x and the binding occurrence of x. This 
variable numbering scheme is often used for implementing functional languages. 
Notice how in this scheme a variable may be denoted differently depending upon 
its position in the term (i.e. depending upon how far it lies from its binding occur- 
rence). The rules of figure[9]implement this mechanism. 

Thus there are two ways to use the A-calculus to define a linear map in a vecto- 
rial space of countable dimension over a countable field. The first is to interprete 
the A-terms as functions mapping vectors to vectors and in this case we need extra 
constraints to enforce linearity. The second is to interprete the A-terms as func- 
tions mapping base vectors to base vectors and to extend it to the full space with 
extra computation rules. This solution is advantageous as it requires no restriction 
on the A-terms. 

Notice that our language allows a restricted form of higher-order programming 
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Figure 9: Explicit substitution of de Bruun indices 



(t + u) bof s — ► 


(t bof s) 


+ (ubof s) 


(r.u) bof s — ► 


(r bof s) 


.(u bof s) 


(t • u) bof s — > 


(t bof s) 


• (u bof s) 


(r + s) bof s — ► 


(r bof s) 


+ (sbofs) 


(r x s) bof s — > 


(r bof s) 


x (sbofs) 


r bof s — ► 


r 




(t <g> u) bof s — > 


(t bof s) 


® (u bof s) 


(t > u) bof s — ► 


(t bof s) 


> (u bof s) 


(t * u) bof s — ► 


(t bof s) 


* (u bof s) 


L(t)bofs — > 


L(tbof ft(s)) 



bof s — > 
false bof s ► false 
true bof s > true 



var(O) bof subst(v) — ► v 
var(S(p)) bof subst(v) — > var(p) 
var(O) bof fr (s) — > var(O) 
var(S(p)) bof fr (s) — > (var(p) bof s) bof | 
var(p) bof | — ► var(5(p)) 



where a function F defined by a A-term can take another function g as its argu- 
ment, provided the function g is expressed with the matching construct, but not if 
it is expressed as a A-term also. The extention of this language to full higher-order 
programming is left for future work. 

6 An example 

The Deutch-Jozsa algorithm can be defined in our language as follows 

DJ = Xx ({Cross * H * H ) * (x * ({Cross * H * H) * (false ® true)))) 
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Figure 10: Summary 
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where 

Cross = Xx Ay ( 

(false <g> false) E> ((x * false) ® (y * false)) 
+ (false ® true) t> ((x * false) £g> (y * true)) 
+ (true (g) false) t> ((x * true) <g> (y * false)) 
+ (true (X> true) l> ((x * true) (y * true))) 

i.e. using de Bruijn indices 

DJ = L((Cross * H * H)* (var(O) * ((Cross * H * H) * (false <g) true)))) 
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where 



Cross = L(L( 

(false <g> false) > ((var(S(0)) * false) ® (var(O) * false)) 
+ (false <g> true) > ((var(S(0)) * false) ® (var(O) * true)) 
+ (true <g> false) t> ((var(5(0)) * true) <g> (var(O) * false)) 
+ (true <g> true) > ((var(5(0)) * true) ® (var(O) * true)))) 



Then it can be checked that the term DJ* CNOT reduces to true <g> true whose 
first component is indeed the exclusive disjunction of not (true) and not (false). 



7 Discussion 
7.1 Unitarity 

In its simplest formulation quantum theory only allows unitary evolutions, i.e. 
vectors evolve in time according to square matrices U verifying Wll = I. In this 
framework it is impossible to delete, say, a qubit: 

|0) E ^ E |0) 
|1) b ^ E |0) 
H0)+/3|l)) b ^ E |0), 

the evolution is not injective and therefore not unitary. Von Neumann's projective 
measurements help us only partially: if the vector is measured in the canonical 
basis, and when this measurement yields outcome "0", then the qubit undergoes 
the above exact dynamics. But this will only occur with probability \a\ 2 . 
The ERASE operation is however perfectly physical, as one can always ignore a 
qubit and focus upon another, taken to be in state |0). Moreover the process needs 
not be probabilistic. There are two well-established formulations of quantum the- 
ory which cater for this possibility: 

• The generalized measurement formalism unifies quantum evolutions and 
quantum measurements as one single object. Mathematically a generalized 
measurement is given by a set of matrices {M m } verifying 



J2MlM m «Id. (2) 
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A vector v will then evolve in time according to the matrix M m with prob- 
ability p m = \M m v\ 2 — in which case we shall say that outcome "m" has 
occurred. As an example the following generalized measurement performs 
the ERASE operation: 



For a more detailed presentation of these concepts the reader is referred to 
lfT8l - page 84. 

• The density matrix formalism represents quantum states as positive matrices 
instead of vectors. These evolve in time according to Completely Positive- 
preserving maps, i.e. operations of the form 

p ' — > ^ MmpM^ with probability p m = Tr(M my oM^) 



and verifying Equation (0. This framework is traditional when dealing 
with open quantum systems, and therefore well appropriate as one discards 
a qubit. 

In classical functional languages terms such as Xy.(Xx.x), with (Xy.(Xx.x)t) -^-> 
Xx.x, are commonly used. Boolean values and branching, for instance, are en- 
coded in such manners. Although convenient, non-injective functions are not 
absolutely necessary, for reversible computation can be both universal and effi- 
cient |5|. Whether a quantum functional language should allow erasure or not 
must therefore depend upon which of the three above mentioned formulation of 
quantum theory gets chosen. 

If we adopt the simplest formulation of quantum theory, our language must 
be restricted to operation which are a not only linear but also unitary. At first 
sight this seems feasible by imposing the standard U'U = I condition upon the 
matching constructs of section|4] and the relevance condition upon the A-terms of 
|5](i.e. for all term Xx.t the variable x must occur at least once in t). This remains 
a subject for future work. 

7.2 Linearity? 

Linear logic appears in 1131 as a mean to express and prove properties of dy- 
namical systems where the consumption of resources is important. The standard 
example (price updated) is A ="I have 6€", B ="I have a paquet of Gauloises", 
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and the statement A — o B to express the possibility of using up A to obtain B. 
With ® now expressing a conjunction, it is clear one cannot have A —° (B ® 5), 
since this would mean buying two paquets for the price of one. Neither can we 
have A® A —o B: we must get something for our money — at worse the feeling of 
getting cheated. Therefore the rules governing symbols — o, ® differ from those of 
classical logic for =>•, A. Unless there is an abundance of resources (denoted by the 
exclamation mark "!"), in which is case they coincide again: IA — o (B(g>B®. . .). 
Whilst considering this point the father of Linear logic has the following thought 
1141 : "Classical logic appears to be the logic of macro -actions, as opposed to 
linear logic which would be a logic of micro -actions. The unusual character of 
linear logic may therefore be considered similar to the strange character of micro- 
mechanics, i.e. quantum mechanics." . 

Specifications expressed in linear logics can be seen as types for programs ex- 
pressed in linear A-calculus. In the linear A-calculus one distinguishes linear re- 
sources, which may not be copied nor discarded, from nonlinear resources, which 
are denoted by the exclamation mark "!" and whose fate is not subjected to par- 
ticular restrictions. Van Tonder's quantum A-calculus (X q ) is founded upon these 
ideas. As we have mentioned in 15. II he uses this well-established framework in 
order to distinguish quantum resources (treated as linear) from classical resources 
(treated as nonlinear): 

(X q ) t ::= x | Xx.t \ (t t)\c\\t\ Xlx.t 
c::=0\l\H\CNOT\P 

(together with the well-formedness rules of the classical linear 
A calculus) 

(%. q ) (Xx.t s) t[s/x] 
H — >0 + l 
H 1 — >0-l 



Here the well-formedness conditions of the classical linear calculus (which pre- 
vents linear terms from being discarded), together with !-suspension (which stops 
quantum terms from being treated as nonlinear) maintain unitarity throughout the 
reductions. 

The connection between the linear A-calculus and quantum functional lan- 
guages is striking. It comes at a price however: the A g -calculus remains hetero- 
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geneous, i.e. a juxtaposition of quantum resources (linear resources) and classical 
resources (nonlinear resources). In some sense it is twice linear, both in the sense 
of linear A-calculus and linear algebra, and thus in some sense overrestricted. In 
particular it forbids both the cloning of quantum data (which needs be done) and 
the copying of quantum data (which needs not be done). As a consequence its 
control flow remains inherently based upon classical resources. 

The linear-algebraic X-calculus constructed in this paper is homogeneous, i.e. 
it does not draw a line between quantum resources and classical resources (the lat- 
ter are merely thought of as basis states of the former). Moreover it exhibits only 
one notion of linearity, which is that of linear algebra. Thus cloning remains dis- 
allowed (just by the semantics) but not copy. Control flow is still provided as a 
consequence. These results seem to open the way to a linear algebraic interpre- 
tation of linear logic, in the spirit of Girard's Geometry of interaction, although 
much work remains ahead in order to strengthen this connection. 
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